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Gamma-ray bursts (GRBs) appear as the brightest transient phenom- 
ena in the Universe. The nature of the central engine in GRBs is a miss- 
ing hnk in the theory of fireballs to their stellar mass progenitors. Here it 
is shown that rotating black holes produce electron-positron outflow when 
brought into contact with a strong magnetic field. The outflow is pro- 
duced by a coupling of the spin of the black hole to the orbit of the par- 
ticles. For a nearly extreme Kerr black hole, particle outflow from an ini- 
tial state of electrostatic equilibrium has a normalized isotropic emission of 
~ 5 X 1Q'^^{B / Bc)'^{M /IMq)"^ s\v? 6 erg/s, where B is the external magnetic 
field strength, Be = 4.4 x 10^^ G, and M is the mass of the black hole. This 
initial outflow has a half-opening angle 9 > y/B^JSB. A connection with 
fireballs in 7-ray bursts is given. 
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Gamma-ray bursts (GRBs) are now believed to be of cosmological origin, as inferred 
from an istropic distribution in the sky 1^^, redshifts of order unity when detected P, |27|j2^ 
and < V jVmax >= 0.334±0.008 distinctly less than the Euclidean value 1/2 (for long bursts) 
| 38| . The observed emission is well described by the internal-external shock fireball model 
36| , |37|j33| , |3^j26| , pJ^ . Their engines are compact, as indicated by the ratio a = GAE/c^6t, 



where AE is the energy released, G is Newton's constant, c is the velocity of light and 
St ~ 1ms is the observed time-scale of variability. The typical values a = 10~^ — 10~^ 
for GRBs is extremely large compared to other burst-type phenomena, such as supernovae 
including the GRB/SN 1998bw event [jl4| and outbursts in accreting systems such as GRS 
1915+105 |p5| , |23[| , and is close to the maximal value of unity. 

The cosmological distances of GRBs indicate a typical energy release of the order of 
IMq, while the value of a suggests an association with a black hole of a few solar masses. 
Black holes are a natural outcome of the evolution of binaries of young massive stars. Their 
evolution branches out over a diverse spectrum of intermediate states, such as accreting 
neutron stars or black holes with a red giant companion. These branches collectively produce 



either supernovae, their "failed" counterparts |^0[ or hypernovae ||3T|J^, or coalescence of 
binaries made of neutron stars, black holes, white dwarfs or compact He-cores (see, e.g. 
||12||). The outcome of these events is probably associated with the formation of a new black 
hole, and most likely through an intermediate black hole/torus or black hole/disk state 
3^,|SD|. The energetic output from these configurations derives from tapping a fraction of 



the gravitational binding energy of the surrounding matter upon accretion onto the black 
hole, or tapping of the rotational energy of the black hole, which can reach about one 
third of its total rest mass. The associated efficiencies and beaming ultimately determine 



the apparent energies. The GRB990123 event |13| is important in this respect, in that its 
apparent energy was in excess of IMq assuming isotropic emission. 

This Letter focuses on the nature of the central engine. A theory is described for electron- 
positron pair-creation powered by a rapidly spinning black hole when brought into contact 
with a strong magnetic field. The magnetic field is supplied by the surrounding matter as 
in forementioned black hole/torus or disk systems. A rapidly spinning black hole couples 
to the surrounding matter by Maxwell stresses p5| , ^ . This coupling suppresses accretion, 
and in the case of a torus formed from the break-up of a neutron star, is expected to be 



intermittent on the time-scales of 0.15-1. 5s ||4^. In this case, the magnetic field is the 
remnant field of its progenitor, and if the torus remains at nuclear density, it may reach 
magnetic field strengths of up to 10^ ''G by linear amplification [^. The theory is described 
by perturbative calculations about a Wald field in electrostatic equilibrium. Particle creation 
towards electrostatic equilibrium has been considered in a previous analysis [|TH|. Here we 
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pay particular attention to the coupling of the black hole spin to the orbit of the charged 
particles, and find it is important in a continuation of the pair creation process. 

Pair-creation can be calculated from the evolution of wave-fronts in curved spacetime, 
which is well-defined between asymptotically fiat in- or out-vacua. By this device, any in- 
equivalence between them becomes apparent, and generally gives rise to particle production 
I^JTj]. It is perhaps best known from the Schwinger process |2S|,|TD|,^, and in dynamical 
spacetimes in cosmological scenarios Such particle production process is driven primar- 
ily by the jump in the zero-energy levels of the asymptotic vacua, and to a lesser degree 
depends on the nature of the transition between them. The energy spectrum of the particles 
is ordinarily nonthermal, with the notable exception of the thermal spectrum in Hawking 



radiation from a horizon surface formed in gravitational collapse to a black hole [jT9[. There 
are natural choices of the asymptotic vacua in asymptotically fiat Minkowski spacetimes, 
where a time-like Killing vector can be used to select a preferred set of observers. This leaves 
the in- and out-vacua determined up to Lorentz transformations on the observers and gauge 
transformations on the wave-function of interest. These ambiguities can be circumvented by 
making reference to Hilbert spaces on null trajectories - the past and future null infinities 
in Hawking's proposal - and by working with gauge- covariant frequencies. The latter 



received some mention in Hawking's original treatise [0, and is briefly as follows. 

Hawking radiation derives from tracing wave-fronts from J+ to J~, past any potential 
barrier and through the collapsing matter, with subsequent Bogolubov projections on the 
Hilbert space of radiative states on J~. This procedure assumes gauge covariance, by tracing 
wave-fronts associated with gauge- covariant frequencies in the presence of a background 
vector potential A^. The generalization to a rotating black hole obtains by taking these 
frequencies relative to real, zero-angular momentum observers (c/. [^,^), whose world- 
lines are orthogonal to the azimuthal Killing vector as given by (,°'da = dt — {gt:j)l g4><i))d^. 
Then ^'^ ~ dt at infinity and ^""da assumes corotation upon approaching the horizon, where 
Qab denotes the Kerr metric. This obtains consistent particle-antiparticle conjugation by 
complex conjugation among all observers, except for the interpretation of a particle or an 
antiparticle. Consequently, Hawking emission from the horizon of a rotating black hole gives 
rise to a flux to infinity 

cPn _ I r 

dujdt ~ 2^e2^('^-^^)/« + 1' ^' 
for a particle of energy u at infinity. Here, k = 1/4M and Qh are the surface gravity 
and angular velocity of the black hole of mass M, F is the relevant absorption factor. The 
Fermi-level Vp derives from the (normalized) gauge-covariant frequency as observed by a 
zero-angular momentum observer close to the horizon, namely, uj — Vp = ojzamo + = 
uj — uQh + for a particle of charge — e and azimuthal quantum number u, where V is the 



3 



potential of the horizon relative to infinity. The results for antiparticles (as seen at infinity) 
follow with a change of sign in the charge, which may be seen to be equivalent to the usual 
transformation rule uj —uj and v ^ —v. 

In case of = Hawking radiation is symmetric under particle-antiparticle conjugation, 
whereby Schwarzschild or Kerr black holes in vacuum show equal emission in particles and 
antiparticles. For a Schwarzschild black hole, then, the resulting luminosity of (|Tp is thermal 
with Hawking temperature T ~ 10~'^(Mq/M)K, which is negligible for black holes of astro- 
physical size |]32| , |4l|| . The charged case forms an interesting exception, however, where the 



Fermi-level —eV gives rise to spontaneous emission by which the black hole equilibrates on a 
dynamical time-scale [|T5i^ |6|. In contrast, the Fermi-level vVLh of a rotating black hole act- 
ing on neutrinos is extremely inefficient in producing spontaneous emission at infinity 



This is due to an exponential cut-off due to a surrounding angular momentum barrier, which 
acts universally on neutrinos independent of the sign of their orbital angular momentum. 
This illustrates that (|l|) should be viewed with two different processes in mind: nonthermal 
spontaneous emission in response to a non-zero Fermi-level, and thermal radiation beyond 
i- 

Upon exposing a rotating black hole to an external magnetic field this radiation pic- 
ture is expected to change, particularly in regards to Vp and the absorption coefficient F. 
The radiative states are now characterized by conservation of magnetic fiux rather than 
conservation of particle angular momentum, which has some interesting consequences. 

Consider a black hole of mass M and specific angular momentum a, which is brought 
into contact with an external magnetic field B aligned with its axis of rotation. Its lowest 
electrostatic energy state is reached upon accretion of a Wald charge |^7| q = 2BJ, where 



J = aM is the angular momentum of the black hole |^8|. Hereby the source- free Wald 
solution of the vector potential, 

Aa = IbK -(t^- aB) ria, (2) 



2 V2M 
reduces to 

Aa = ^Bka, (3) 

where k"- is the azimuthal Killing vector and rj"- is the asymptotically time-like Killing vector. 
In electrostatic equilibrium, then, ^"A^ = 0. The charge q can be understood by noting that 
it restores the fiux through a hemisphere of the horizon of a rapidly spinning black hole . 
On dimensional grounds, it contributes a fiux $' = kqMVLu = 2A;i?M^ sin^(A/2), where k 
is a dimensionless constant of proportionality and Qh = tan(A/2)/2M using sin A = a/M 
1 ^ . With k = Air, $' plus the fiux AirBM"^ cos A through a charge- free horizon recovers the 
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Schwarzschild value $ = AnBM'^, as given by the exact expression (Q). It should be noted 
that the equilibrium value q = 2BJ obtains the uncharged Kerr metric upon neglecting 
gravitational contributions of the stress-energy tensor of the electromagnetic field [^]. The 
state of electrostatic equilibrium of a black hole in vacuum is therefore described by 
corresponding to a maximal horizon flux 2ttA^ P7|,|9|. 

In an axisymmetric magnetic field B parallel to the axis of rotation, the wave-functions 
of charged particles can be expanded locally in coordinates (p, 0, s,t) as e~*'^*e*'^'^e*^'"^ ■?/'(/)), 
where s denotes arclength along the magnetic field. Comparison with the theory of plane- 
wave solutions gives a localization on the u-th flux surface at which gV? = J2u/eB 



with Landau levels Ena = {"ml + pi + \eB\{2n + 1 — a)}^/^, where rUe is the electron mass 
and a = ±1 refers to spin orientation along B. These states enclose a flux A^p = \Bk'^ = 
v/e. Note that the latter have effective cross sections Sj, = 2'nj\eB\. The gauge-covariant 
frequency of the Landau states near the horizon follows from 

- i'^ii-^da + eAa)^ = {yj- uQh)^. (4) 

The jump 

Vf = i-n^-'da + eA^)]gij = u^H (5) 

between the horizon and infinity defines the Fermi-level of the particles at the horizon. In 
contrast, the Wald field about an uncharged black hole has Vp = v^h — eaBo, which shows 
that it is out of electrostatic equilibrium. Note that the canonical angular momentum of 
the Landau states vanishes: k""n'aip = {i^^d^ — eA^)il) = 0. The Fermi-level combines the 
spin coupling of the black hole to the vector potential Aa and the particle wave-function 
ip. The equilibrium state in the sense of dtq ~ 0, or at most q/dtq ~ a/dta, derives from 
this complete Vp- For this reason, we shall study the state of electrostatic equilibrium as an 
initial condition, to infer aspects of the late time evolution. 

The strength of the spin-orbit coupling which drives a Schwinger-type process on the 
surfaces of constant flux may be compared with the spin coupling to the vector potential 
Aa- The latter can be expressed in terms of the EMF^ over a loop which closes at infinity 
and extends over the axis of rotation, the horizon and the z/— th flux surface with flux 
Thus, we have EMFy = VLh'^,j/2ti PP3[1, which gives rise to the new identity 



eEMF^ = vVLh. (6) 

It should be mentioned that continues to hold away from electrostatic equilibrium (i.e.: 
q 7^ 2BJ), since C,"'Aa = and hence u — eA^ = on the horizon. Since the latter is a 
conserved quantity, it, in fact, continues to hold everywhere in the Wald-field approximation. 
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In the assumed electrostatic equilibrium state, ^""Aa = 0, and the generalization of (H) 
to points (s, z/) away from the horizon is 

[-ar'da + eAa)]^''^^ = -ly^is, u) = -\eBg,4s, u) = -eA.is, u) (7) 



for particles of charge — e. Thus, (|^) locahzes (|^) by expressing the coupling of the black 
hole spin to the wave-functions in terms of the electrostatic potential V = At in Boyer- 
Linquist coordinates. Note that the real zero-angular momentum observers detect no electric 
potential. 

The particle outflow derives from the distribution function (|l]) by calculation of the trans- 
mission coefficient through a barrier in the so-called level-crossing picture @]. The WKB 
approximation (e.g.: as derived by zero-angular momentum observers) gives the inhomoge- 
neous dispersion relation 

(uj-VfY = ml + \eB\{2n + l-a)+pl (8) 

where Vp = Vf{s, v) is the s— dependent Fermi-level on the v-ih flux surface. The classical 
limit of is illustrative, noting that the energy e of the particle is always the same relative 
to the local ZAMOs that it passes. Indeed, since w°'{maa — eAa) is conserved when w"" is a 



Killing vector [^, rj°'{mUa — eAa) = vr^ and k'^{mUa — eAa) = are constants of motion, 
where is the four- velocity of the guiding center of the particle, and tti = Ena, tt^ = in 
a Landau state. With ^""Aa = 0, e = —^°"mua = —^""{mua — eAa) = —rj"-{mua — eAa) = vr^. 
This conservation law circumvents discussions on the role of E ■ B (generally nonzero in a 
Wald field). The energy of the particle relative to infinity is uj. This relates to the energy 
e as measured by the ZAMOs following a shift Vf{s, v) due to their angular velocity. Thus 
(P) pertains to observations in ZAMO frames, but is expressed in terms of the energy at 
infinity uj. It follows that particle/antiparticle pair creation (as in pair creation of neutrinos 
is set by 



?7 = \dVF/ds\ 



dr ( ]^eBgt4> 



\dAeM)\ = eBaM—^^—^^ sin^ ^, (9) 



using ds ^ dr- Radiation states at infinity are separated from those near the horizon by a 
barrier where < about Vf{sq) = uj. The WKB approximation gives the transmission 
coefficient 

Since the Wald field B is approximately uniform, any additional magnetic mirror effects can 
be neglected. Also, r] < \eB{M / a)iaxi^ 6 < \eB and \eB\{2n + I - a) / -q > A(2n + 1 - a), 
so that T is dominated by n = and a = 1. 
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By the pair production rate by the forcing 77 in (D can be derived from the analogous 
results for the pair production rate produced by an electric field E along B. The results 
from the latter |^,^ imply a production rate of particles given by 

47r2 y tanh(7re5/77)^ ^ 2ti J + cos^ 6 ^ ' 

Here l/r] ~ {eBaM sin^ d)-^{9>a^ + 12(r - VSa cos^)^) about r = -\/3flcos^. For a rapidly 
spinning black hole, -y/Sacos^ is outside the horizon in the small angle approximation, 
whereby after r-integration of (|l^) we are left with 

^ ^ eWM /g-8Wsin^^in4^rf^ ^ ( ^\\-V^e-^-c/e-e' (12) 

asymptotically as '^tic/O'^ » 1. Here c = mla/eBM , Nh = m^M'^ is a characteristic 
number of particles on the horizon, and 9 is the half-opening angle of the outflow. The 
right hand-side of (p!2D forms a lower limit in case of Sttc/O < 1. When a ~ M, Nh/c 
is characteristic for the total number of flux surfaces which penetrate the horizon and 
c ~ Bc/B, where Be = 4.4 x lO^^G is the field strength which sets the first Landau level at 
the rest mass energy. By (^) and (^^, a similar calculation obtains for the luminosity in 
particles Lp normalized to isotropic emission the estimate 

~ ^ ~ VSeBMN. (13) 

This analysis shows that the black hole departs from electrostatic equilibrium by outflow 
of e~ (with the sign convention BQ^ > 0) towards q > 2BJ. The calculations therefore 
pertain to the initial jet formation. Upon accumulating q > 2BJ, the jet evolves towards 
an inner jet of outflow produced by Vp which is dominated by a positive electrostatic 
potential near the polar caps, surrounded by e~ outflow produced by Vp which is dominated 
by the spin-orbit coupling as in the initial jet. A full calculation of the resulting equilibrium 
outflow (in the sense as described before) falls outside the present scope. Nevertheless, 
it is expected that the luminosity (^3|) remains characteristic for the combined particle- 
antiparticle outflow in the evolved jet. 

The outflow ( p!3D as follows from (0) saturates when the back reaction of the outflow 
is taken into account. This is non-dissipative in the form of an attenuation of the mag- 
netic field by the collective contribution of the charged particles to azimuthal currents, and 
dissipative due to a finite surface conductivity of An (see |^) in current closure over the 
horizon. The four-current of the charged particles, and the associated a Poynting flux, forms 
a perturbation away from the source-free Wald field in our derivation. While these two types 
of back reactions are similar in order of magnitude, in angular dependence the dissipative 
back reaction is dominant. It follows that 
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AireN < vVLh 



(14) 



for the outer, spin-orbit driven e~ outflow to proceed, up to a logarithmic factor of order 
In {tt/26), where u is considered at the half-opening angle 6 of the outflow. Consistency with 
(p!^) gives for the minimum opening angle Oq for the outflow in e~ 



up to logarithmic corrections. For 6 > Oq, the outflow is effectively set by the saturation 
limit (0), whereby the particle luminosity ([T3|) becomes 

where ^ > is the half-opening angle of the outflow. This calculation applies formally to 
the initial jet. A full calculation of the evolved jet, which consists of combined e^-outflow 
saturated against dissipative losses in the horizon (in the sense as described above) falls 
outside the present scope. Nevertheless, it is expected that the luminosity (|16D remains 
characteristic for particle-antiparticle outflow in the evolved jet, whose opening angle will 
be bounded below by the initial value (p^Sf). 

A connection to fireballs p6|j37|JTl| in the theory of 7-ray bursts is at hand, as (|T6|) 



represents a GRB type of luminosity for B ~ lO^^G, whose energetics are consistent with 
the GRB990123 event, with ultra- relativistic electron-positron outflow along open field-lines 
(supported by the horizon charge q) in agreement with the input of current fireball models. 

The theory of fireballs accounts for the high luminosity in nonthermal emission either by 
way of baryon loading or intermittency, to circumvent the thermal emission in the preceding 



steady-state models of electron-positron fireballs |T^,[18[. A small amount of baryonic matter 



enables efficient conversion of the energy in the fireball into kinetic energy, with subsequent 



shocked emission in interaction with the interstellar medium or wind [|39| , p6| . Some contam- 
ination should be expected, for example, when the interstellar medium or a wind from the 
surrounding torus is entrained, and, in case of hypernovae, by interaction with the hydrogen 



envelope. On the other hand, intermittency at the source (see e.g. [^,^,0) will give rise 



to unsteadiness in the flow as discussed in the compact fireball model of Eichler & Levinson 
||Tl| with angular variations (both in Oq and in orientation) and internal shocks even when 
baryon free. Shocks may also result from interactions with coUimating baryonic material, 
perhaps subject to radiative viscosity, whereby a nonthermal component in the emission is 
produced from a broad range of radii [|1T| . These considerations imply an emission spectrum 



substantially different from that in forementioned steady-state models. 
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